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We call projectile fragmentation of neutron halo nuclei the elastic breakup (diffraction) 
reaction, when the observable studied is the neutron-core relative energy spectrum. This 
observable has been measured in relation to the Coulomb breakup on heavy target and 
recently also on light targets. Such data enlighten the effect of the neutron final state 
interaction with the core of origin. Projectile fragmentation is studied here by a time 
dependent model for the excitation of a nucleon from a bound state to a continuum 
resonant state in a neutron-core complex potential which acts as a final state interaction. 
The final state is described by an optical model S-matrix so that both resonant and non 
resonant states of any continuum energy can be studied as well as deeply bound initial 
states. It turns out that due to the coupling between the initial and final states, the 
neutron-core free particle phase shifts are modified, in the exit channel, by an additional 
phase. Some typical numerical calculations for the relevant observables are presented 
and compared to experimental data. It is suggest that the excitation energy spectra 
of an unbound nucleus might reflect the structure of the parent nucleus from whose 
fragmentation they are obtained. 

1. Introduction 

All theoretical methods used so far to describe breakup rely on a basic approximation 
to describe the collision with only the three-body variables of nucleon coordinate, pro- 
jectile coordinate, and target coordinate. Thus the dynamics is controlled by the three 
potentials describing nucleon-core, nucleon-target, and core-target interactions. In most 
cases the projectile-target relative motion is treated semiclassically by using a trajectory 
of the center of the projectile relative to the center of the target R (t) = b c + vtz with 
constant velocity v in the z direction and impact parameter b c in the xy plane. This 
approximation makes our formalism applicable for incident energies above the Coulomb 
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barrier. Along this trajectory the amplitude for a transition from a nucleon state ipi bound 
in the projectile, to a final continuum state is given by [[TJ [2] 



where V is the interaction responsible for the transition which will be specified in the fol- 
lowing. The probabilities for different processes can be represented in terms of the ampli- 
tude as dP/d^ = J2 IAfi| 2( K£ — £/) where £ can be momentum, energy or any other variable 
for which a differential cross section is measured. Direct one-particle re-arrangment reac- 
tions of the peripheral type in presence of strong core-target absorption can be described 
by an equation like [ [2j [3j HI [5] 



(see Eq. (2.3) of [[3]) and C 2 S is the spectroscopic factor for the initial single particle 
state. The core survival probability is defined in terms of a S-matrix function of the core- 
target distance of closest approach b c . A simple parameterisation is P c t(p c ) = \S ct \ 2 = 
e (- in2exp{(R a -b c )/a]) _ j^. t a k es [ n ^ acc0 unt the peripheral nature of the reaction and natu- 
rally excludes the possibility of large overlaps between projectile and target. The strong 
absorption radius R s ~ 1A(A^ 3 + A^ 3 ) fm is defined as the distance of closest approach 
for a trajectory that is 50% absorbed from the elastic channel and a=0.6 fm is a diffuse- 
ness parameter. The values of R s thus obtained agree within a few percent with those of 
the Kox parameterization [[6]. 

2. Projectile Fragmentation 

Let us call projectile fragmentation the elastic breakup (diffraction dissociation) reac- 
tion, when the observable studied is the neutron-core relative energy spectrum. This kind 
of observable has been widely measured in relation to the Coulomb breakup on heavy 
target. Results on light targets have also been presented [[?]. These data enlighten the 
effect of the neutron final state interaction with the core of origin, while observables like 
the core energy or momentum distributions enlighten the effect of the neutron final state 
interaction with the target. 

Projectile fragmentation has also been used to study two neutron halo projectiles [[9]-[ 
[21]. In this case it has been suggested that the reaction might proceed by the simultaneous 
emission of the two neutrons or by successive emissions [ [9] . The successive emission can 
be due to a mechanism in which one neutron is stripped by the interaction with the target, 
as in the one-neutron fragmentation case, while the other is left behind, for example in a 
resonance state, which then decays. This mechanism has been described by the sudden 
approximation [ [10] under the hypothesis that while the first neutron is stripped, the 
second neutron is emitted at large impact parameters with no final state interaction with 
the target. The emission can be expected sequential if the two neutrons are not strongly 
correlated. 

If the two neutrons are strongly correlated they will preferentially be emitted simul- 
taneously. If the neutron which is not detected is stripped while the other suffers an 
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elastic scattering on the target, then in both cases to first order in the interaction the 
neutron ends-up in a plane wave final state [ 1\ . It can then re-interact with the core 
which, for example, is going to be 10 Be in the case of the one-neutron halo projectile 11 Be, 
while it will be 12 Be in the case of the projectile fragmentation of 14 Be, since 13 Be is not 
bound. Experiments with a 14 B projectile [ [11] have also been performed, in which the 
n- 12 Be relative energy spectra have been reconstructed by coincidence measurements. In 
such a nucleus the valence neutron is weakly bound, while the valence proton is strongly 
bound. Thus the neutron will probably be emitted in the first step and then re-scattered 
by the core minus one proton nucleus. The projectile-target distances at which this kind 
of mechanism would be relevant are probably not so large to neglect the effect of the 
neutron-target interaction. 

2.1. Inelastic excitation to the continuum. 

To first order the inelastic-like excitations can be described again by the time dependent 
perturbation amplitude Eq.([T]) [HI [2]- In this section also, the potential ^(r, R(t)), which 
is the interaction responsible for the neutron transition, moves past on a constant velocity 
path as described in the previous sections. The radial part 0«(r) of the single particle 
initial state wave function ipi(r,t) is calculated in a potential Vws{ r ) which is fixed in 
space. The coordinate system and other details of the calculations can be found in Ref. [ 
[T2] . In the special case of exotic nuclei the traditional approach to inelastic excitations 
needs to be modified. For example the final state can be eigenstate of a potential V\ 
modified with respect to Vws because some other particle is emitted during the reaction 
process as discussed in the introduction. The final state interaction might also have an 
imaginary part which would take into account the coupling between a continuum state 
and an excited core. The first order time dependent perturbation amplitude then reads 

A fi = — dxdydz (f)* f (x,y,z)(j) i (x,y,z)e iqz V(x -b c ,y,q), (3) 

where V(x — b c , y, q) = dzV(x — b c , y, z)e tqz , and we changed variables and put z' = 
z — vt or t = (z — z')/v, q = Ef — Ei/hv. Here £/ is the neutron-core relative energy in 
the final state. 

The role of the target represented by V is just to perturb the initial bound state wave 
function and to allow the transition to the continuum by transferring some momentum to 
the neutron. Then it is enough to choose a simplified form of the interaction, such as a 
delta- function potential V(r) = V2S(x)5(y)5(z). The value of the strength v 2 = [MeV fm 3 ] 
used in the calculation is taken equal to the volume integral of the appropriate neutron- 
target interaction. It is clear that while in the sudden approach the initial and final state 
overlap is taken in the whole coordinate space, irrespective of the target and of the beam 
velocity, here the overlap of the initial and final wave functions depends on the core-target 
impact parameter. The neutron is emitted preferentially on the reaction plane and the 
z-component, being along the relative velocity axis is boosted by a momentum q. 

Due of the strong core absorption discussed in Sec.l these calculations are performed 
using the asymptotic form of the initial and final state wave functions. Introducing the 
quantization condition[[2] the probability spectrum reads 

dPin _ 2 V I „2 m \ y I -I n |2|r |2 ( a\ 

de f ~ 7rh 2 v 2 1 h 2 k2k + 1 mi>m/l m - mfl 1 mi ' mfl ■ { ' 
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The generalization including spin is given in Appendix B of Ref . [ [T2] and \I. 



/^d^e^T^fiTr)^,^^, 0)^1^ (Ar)^ /)m/ (0,0) ". The quantity S = e 2i( - 6+ ^ is an 
off-the-energy-shell S-matrix representing the final state interaction of the neutron with 
the projectile core. It depends on a phase which is the sum of S, the free particle n-core 
phase shift, plus v the phase of the matrix element |/|. 
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3. Applications 

3.1. The reaction n Be -> n+ 10 Be 





Figure 1. n- 10 Be relative energy spectrum, including Coulomb and nuclear breakup for the reaction 
11 Be+ 12 C — ► n+ 10 Be+X at 69 A. MeV. Only the contributions from an s initial state with spectroscopic 
factor C 2 S= 0.84 are calculated. The triangles are the total calculated result after convolution with the 
experimental resolution function. The dots are the experimental points from [[7j. 



As a test of our model we calculate the relative energy spectrum n+ 10 Be obtained by 
the authors of Ref. [UJ in the breakup reaction of n Be on 12 C at 70 A. MeV. The structure 
of n Be is well known: the valence neutron is bound by 0.503 MeV; the wave function is 
mainly a 2s state with a spectroscopic factor around 0.8 and there is also a small ds/2 
component. The main d.5/2 strength is in the continuum centered around 1.25. We have 
calculated the initial wave function for the s-state in a simple Woods-Saxon potential with 
strength fitted to the experimental separation energy and whose parameters are: ro=1.25 
fm, a=0.8 fm. As possible final states we have considered only the s, p and d partial waves 
calculated in the /-dependent potentials of [H2]. The delta-function potential strength 
has been chosen as -4057.59 MeV fm 3 . The authors of Ref.[[7J have shown that the effect 
of Coulomb breakup is noticeable in their n+ 10 Be spectrum. We have also included this 
contribution, calculating it according to [Hj. The spectrum of Fig l3.ll is very similar to 
the spectrum obtained in Ref. [ [8] by solving the time-dependent Schrodinger equation 
numerically, expanding the projectile wave function upon a three-dimensional spherical 
mesh. Similarly to the present model, a classical, straight line trajectory for the core- 
target scattering was used in Ref.[|8]. Also our n-core potentials are very close to those 
of Ref. [ [8] and our 5-interaction strength is consistent with the volume integral of their 
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neutron-target interaction. We have then folded the calculated spectrum througth the 
experimental resolution function of Fukuda et al. [ 7], as given in Ref.[[S]. The result is 
shown in Fig |3.1l by the triangles. The full curve is the total spectrum, sum of Coulomb 
and nuclear breakup. Each individual transition, due to the nuclear interaction only, is 
also shown. The dots are the experimental points from [[7j. The kind of discrepancy 
between our calculation and the data in the range 1-2 MeV is very similar to that of the 
calculations in Ref . [ [S] . 

3.2. Structure of 14 Be and 14 B 

Uncertainties in the interpretation of experimental results as compared to structure cal- 
culations were at the origin of our motivations to try to understand whether the neutron- 
12 Be relative energy spectra obtained from fragmentation of 14 Be or 14 B would show dif- 
ferences predictable in a theoretical model. If differences will be found in the experimental 
results with 14 B and 14 Be beams they could be due to an interplay between structure and 
reaction effects. 

The ground state of 14 Be has spin J n = + . In a simple model assuming two neutrons 
added to a 12 Be core in its ground state the wave function is: 

\ lA Be >= [b x {2s 1/2 f + b 2 (l Pl/2 ) 2 + b 3 (ld 5/2 ) 2 ] ® \ 12 Be,0 + > . (5) 

Then the bound neutron can be in a 2s, lpi/2 or ld 5 / 2 state. However, as it has been 
discussed in the previous section, the situation is much more complicated [ [18]-[ [22] 
and in particular the calculations of Ref. [ [23] show that there is a large component 
(2si/2, 1^5/2) ® \ 12 Be,2 + > with the core in its low energy 2 + state which can modify the 
neutron distribution. 




e f (MeV) 



Figure 2. Sum of all transitions from the s initial state for the reaction 14 Be+ 12 C — » n+ 12 Be+X . Ex- 
perimental points from [US]. Dashed line is the folding of the calculated spectrum with the experimental 
resolution curve. 



The ground state of 14 B has spin J w = 2~ . In a model where it is described as a 
neutron-proton pair added to a 12 Be core in its + ground state with the proton in the 
lp3/2 shell, its wave function may be written as: 

\ U B >= [ai(p3 /2j 2s 1/2 ) + a 2 (p 3/2 , d 5/2 )} <g> \ 12 Be, 0+ > . (6) 
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The present experimental information [ 25J on 14 B is that the neutron is in a state combi- 
nation of s and d-components with weights 66% and 30% respectively, while shell model 
calculations show a similar mixture and no component with an excited state of the core. 
There are two possibilities for the reaction mechanism. One is that a proton is knocked 
out in the reaction with the target. The remaining 13 Be would be left in an unbound 
s-state with probability |ai| 2 , in a d 5 / 2 -state with probability |a 2 | 2 . These unbound states 
would decay showing the s-wave threshold and d-wave resonance effects. As mentioned 
in the introduction, the second possibility is that the neutron is knocked out first due to 
its small separation energy and that the proton is stripped from the remaining 13 B. 

To give another example of a possible comparison with available data, we show in Fig. 
13.21 the experimental points from H. Simon et al. [ [16] for the reaction 14 Be+ 12 C — > 
n+ 12 Be+X at 250 A.MeV. The normalization factor of the data to mb/MeV is 0.843. 
The solid line gives the sum of all transitions from the s initial state with £y=-1.85 MeV 
(solid line), renormalized with a factor 2.4. The dashed line is the folding of the calculated 
spectrum with the experimental resolution curve. Therefore the calculation underestimate 
the absolute experimental cross section by a factor of 2. In view of the incertitude in the 
strength of our n-target (^-potential and on the initial state spectroscopic factor which has 
been taken as unit, we can consider our absolute cross sections quite reasonable. A more 
detailed account of these calculations is given in [[T2"]. 

4. Conclusions and Outlook 

The field of Rare Isotopes Studies is very active, growing steadily and rapidly. Some 
recent achievements in the reaction theory for elastic breakup have been presented. From 
the structure point of view, in the search for the dripline position, a very important role 
is played by the study of nuclei unstable by neutron emission. This is one of the most 
important subjects which need to be adressed and further developed in the near future 
and for which some suggestions have been presented. 
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